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A toroidal current distribution has nonvanishing exterior vector potential A, but zero exterior field 
B=VxA=0. This property, together with the absence of fringing fields as in a cylindrical solenoid, 
makes it convenient for studies involving the vector potential in a field-free region, such as the 
Aharonov-Bohm effect, or the effect of A in a Josephson junction. We present an immediate general 
result of magnetostatics, and use it to easily compute A for a torus, to visualize the static vector 
potential for any current distribution, and to show how one can construct a current distribution to 
produce any desired A. When the torus current I varies in time, nonzero quasistatic fields E (t) and 
B(f) are produced (E —tol/r 3 and B ~<u 2 //r 2 ). Radiation is also produced, with the radiation pattern 
of an electric dipole. The torus provides a counterexample to the common erroneous notion that if 
all multipole moments of a current distribution vanish then quasistatic fields and radiation must also 
vanish. We then formulate Maxwell’s equations in a way that obviates the role of gauge 
transformations. This “gauge irrelevant” form clarifies the relation of potentials to current sources, 
isolating the role of the transverse part of A. The general result from magnetostatics is extended to 
time-varying sources, revealing a seldom recognized symmetry of Maxwell’s equations, and 
showing how one can visualize A for an arbitrary time-dependent current source. © 1995 
American Association of Physics Teachers. 

I. INTRODUCTION 

In classical electromagnetic theory, the magnetic vector 
potential A is generally considered only a mathematical aid 
to solving for the field B=VxA. It is the field that is real. 1 

However, in classical physics energy is quite real, and the 
electrostatic potential <f> is usually considered as real as the 
electric field E=-V<£. Since relativity requires r/> and A be 
components of a four vector, one should attribute just as 
much reality to A. Moreover, in quantum mechanics, the 
canonical quantization procedure requires the use of poten- 
tials. A problem would seem to arise in those situations in 
which B=0, but A¥=0, for there should be no classical dif- 
ference, but there may be a QM difference. In spite of these 
observations, the classical belief in the preeminent position 
of the fields held sway until nearly 1960. 

Aharonov and Bohm 2 first pointed out experiments to 
demonstrate the reality and importance of the potentials in 
quantum mechanics. There are observable differences when, 
say, an electron is passed through a region of zero fields but 
nonzero potentials. The differences show up in the phase of 
the wave function, requiring an electron interference experi- 
ment to detect. 

Ever since the reality of A in quantum mechanics was 
emphasized, and various experiments confirmed it, 3-7 inter- 
est has attached to measuring A directly, 8 especially where 
B=0. Patents have been issued for a vector potential mea- 
suring device, 9 and a recent proposal suggests a nondestruc- 
tive photon detector by “passively” measuring only A of the 
photon. 10 Indeed the possibility then arises of a class of elec- 
tromagnetic sensors that would detect the vector potential 
rather than the fields. 

To investigate this phenomenon, it is useful to have a 
ready current distribution that provides a working volume in 
which B=0, but A#0. This is the case for the static A out- 
side a torus on which flows a steady current only in the 
“toroidal” direction, around the thin limb. (But it is not the 
case for time varying currents.) 

Two suggested detectors of A both operate via variants of 
the Aharonov-Bohm effect. In place of a coherent electron 


beam split into two sub-beams enclosing magnetic flux, one 
can use a superconductor with a Josephson junction. The 
superconducting state is a macroscopic coherent wave func- 
tion that plays the role of the electron beam. The Josephson 
junction detects A since the tunneling current depends on the 
phase of the wave function, which shifts when A^O. This is 
the basis of patents by Gelinas. 9,11 

A second detector 10 would employ the conventional split 
electron beam along the surface of a crystal, detecting the 
evanescent A of a light beam undergoing total internal 
reflection. 12 

These sensors work best on a time varying A. Therefore 
for time dependent laboratory current sources, one needs to 
know the fiill environment at the detector, both fields and 
potentials. 

A torus provides a simple source of nonzero A but zero B 
in a workable volume of space. In this note we first obtain a 
very useful expression for the static A outside a torus, and 
for the quasistatic and radiated fields and potential when the 
current varies in time. An Appendix presents the exact fields. 

In the process we make useful observations in magneto- 
statics and electrodynamics concerning the calculation of the 
vector potential, and its relationship to sources in the light of 
gauge invariance. Maxwell’s equations are formulated in a 
way that explicitly incorporates gauge invariance; once in 
that form questions of gauge transformations and gauge in- 
variance never arise. 


II. MAGNETOSTATICS 

The equations of magnetostatics, 

V B=0, VxB=/r 0 J, (2.1) 

are conveniently formulated with 

B=VxA, V -A arbitrary, (2.2) 

so that A obeys 

Vx(VxA) = V(V-A)-V 2 A=u 0 J. ( 2 . 3 ) 
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The arbitrariness in V-A means the gradient of any scalar 
may be added to A without changing B (gauge transforma- 
tion); it is sufficient to compute A in any gauge. Noting that 
in the static case the Lorentz (V-A+ <?</>/<?! =0) and Coulomb 
(V-A=0) gauges are the same, we choose 

V -A=0. (2.4) 

Then A obeys Poisson’s equation V 2 A=-poJ, with the so- 
lution 


. , , Mo f , 3 , J(r ) 

4 77 J | r— r | 

However the basic equations of A, 
V A=0, V xA=B, 


(2.5) 

( 2 . 6 ) 


are the same as obeyed by B, Eqs. (2.1), with A replacing B 
and B replacing p^J. Therefore A is the same function of B 
as B is of yUoJ. We can introduce a vector A by 


A= VxA 


(2.7) 


and choose V-A=0, so that 


( 2 . 8 ) 


F 4rr J r /4 0 |r-r'|' 

Taking the curl of this shows A is given by a Biot-Savart 
law in terms of B, 

a Q f , r' — r B(r') 
AM-VXAM-^J^ x— . (2.9) 


Thus, as an immediate result of magnetostatics, B is the 
source of A just as J is the source of B. This simple obser- 
vation, often overlooked, can be used to visualize A in any 
situation, for A is constructed from B in the same way B is 
constructed from J. One can infer features of A without hav- 
ing to carry out a calculation; for example, the “right-hand 
rule” determines A’s direction from that of B. And A can be 
computed quantitatively in two (identical) steps: first B from 
J, then A from B. 

Magnetostatics admits an endless hierarchy of potentials, 


VxJ= (prescribed) V-J=0, 

V xB=J V B=0, 

V xA=B V -A=0, (2.10) 

V xA=A V -A=0, 


each determined from the next one by the curl operation. B, 
for example, is the “potential” of the “field” J, and A is the 
“field” of the “current” B. This can be very useful in visu- 
alizing and computing the vector potential. 

The hierarchy shows that to obtain the vector potential A x 
of a current distribution J x with field B lf we imagine a sec- 
ond current J 2 equal to B x . Then Aj is equal to the field B 2 of 
J 2 . For any J the solution for B immediately provides the 
solution for an infinite set of problems, and the hierarchy 
may be summarized as 




Fig. 1. Torus geometry definitions. 


This is a general symmetry of magnetostatics. Use will be 
made of these observations in Sec. IX to construct a hierar- 
chy of current distributions based on that of a torus. Equation 
( 2 . 11 ) also allows easy construction of many current distri- 
butions which produce vanishing B but nonvanishing A . In 
Sec. X, it is extended to time-varying sources and fields. 


IH. STATIC VECTOR POTENTIAL OF TORUS 


Direct evaluation of (2.5) for a torus runs into integrals of 
elliptic integrals. While these cannot be avoided for an exact 
solution by quadratures, the observations of Sec. II permit a 
very useful approximate expression. 

The torus lies in the x,y plane, has minor radius a and 
major radius b (Fig. 1). To keep B=0 outside, we require no 
circumferential current in the azimuthal direction about the 
symmetry axis z. The only allowed current is in the direction 
of increasing a, as would be created by a tightly spaced 
toroidal wire coil with an even number of counter-rotating 
layers. If i w is the wire current, and the total number of turns 
is N, then the total current is I=Ni w . 

We employ the usual spherical coordinate system ( r,0,<P ), 
a Cartesian set of axes ( x,y,z ), the usual cylindrical system 
(p,<p,z), and occasionally the internal polar coordinates (s,a). 

The field inside is 


B=B^= 


Vo I 
27 rp ’ 


(3.1) 


and decreases across the interior. The surface current density 
is X=//27rp=//27r(h + acos a) and the current density J is 
KS{s-a)a. The flux in the torus is 



dS-B=7ra 2 


-el-} 

27 Tb 8 \bj’ 


(3.2) 


where dS=sds da, and, with u = s/a, ij=a/b, 

1 _( 1 _£ 2 ) 1 /2 


ri r 7r da 1 

g(£) = 2 duu — — — - = 2 

Jo Jo ir l + «f cos a 


(3.3) 


is the shape factor, l^g^ 2 . g — >1 as a/b—>0 (bicycle tire), 
and g— ► 2 as a/b—> 1 (holeless donut). Torus self inductance 
is L = /IqN 2 (a 2 / 2b) g. 

In the gauge V-A=0, the vector potential is given by the 
curl of ( 2 . 8 ), 


__ Vo . 7 , 

A= V x- — d 3 r' 

477 


y'B(r') 

Mok-r' 


(3.4) 
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Fig. 2. Lines of A. 


the integral being taken over the torus volume. A has no <p 
component. With B given by Eq. (3.1), this equation shows 
A is exactly the same as the magnetic field of a single fat 
wire loop coinciding with the torus with current distribution 
inside the wire proportional to 1/p. When b>a, or when the 
observer is at r>(b + a), this current variation within the 
wire is not important, and A looks like the dipole field of an 
ordinary current loop. Indeed the torus current density J is 
the curl of the current density of an ordinary loop (of wire 
radius a), so (2.11) applies. 

The integral in Eq. (3.4) behaves as 1/r 2 for r>(b+a), 
and there is easily evaluated to obtain 

H o VI no VI 

Af = 4 tt 2ttt^ COS 


proportional to torus volume V=2Tr 2 a 2 b. L ines of A are 
sketched in Fig. 2. The exact solution in Appendix A shows 
that A contains only odd powers of 1/r. The next correction 
to the static potential (3.5) is ~l/r 5 , smaller by a factor 
~(b/r) 2 . 

Outside the torus, A is locally the gradient of a scalar, 
A= — V'P [for the 1/r 3 terms, (3.5), 'P is 

Mo VI 

(3 ' 6) 

Thus A can be locally transformed away with the gauge 
transformation 

A— >A'=A+V'P = 0. (3.7) 

However A cannot be transformed to zero everywhere in 
the doubly connected space outside the torus. Due to (2.2), 
we have around any closed path encircling the limb, 

<j) <fl-A=<l>. (3.8) 

Therefore, on the surface, the average A around the limb is 
(i 0 al 

A ~2TTa~4n'b 8, ^ 

in any gauge. The largest A occurs for a = b, 
A max (mo/47t)2/ . 

In general, since A*$>, but B°c/, A outside can be in- 
creased relative to any leakage B by keeping I constant and 
increasing a (so long as the gap between wires is not also 
increased). 

The equivalence of A of a torus to B of a current loop 
becomes clear by dividing Eq. (3.4) by a, and using Eq. 
(3.1): 


-A =^VxfdV- 7 r. (3.10) 

a 47 t J 2irap|r— r | 

A of a torus is a times the magnetic field of an ordinary loop 
with current density 


^ 27 rap ^ 

and current 


A. Potential on axis 

On the z axis, A is simple when a<b. Then the magnetic 
field of a ring current ( a/2b)I is 13 

fi 0 2vb 2 (al2b)I 

(3 ‘ 13) 

so that A z for a torus with a<b and current I is 

Un 77a 2 W 

(3 ‘ 14) 

A? is largest at z=0, and there is quite comparable to the 
distant components equation (3.5) extrapolated back to r^b, 
but smaller by about a factor ira/b than A on the surface. 


B. The Aharonov-Bohm effect 

Equation (3.8) is the important relation for the Aharonov- 
Bohm effect. The phase shift of an electron moving from 
point P to point Q is 2 -(e/fi)/pdl-A. An electron travers- 
ing a path 1, which passes through the torus, suffers a differ- 
ent phase shift Stp from one traversing path 2, which does 
not, by a gauge invariant amount 


(M5) 

even though fields are zero along both paths. Therefore, the 
interference pattern at Q of a coherent electron source ema- 
nating from P, passing through the two slits with one torus 
arm between them, differs according as the magnetic flux 
between paths 1 and 2 is zero or nonzero, even though no 
force acts on an electron taking either path. 


IV. ALTERNATE CONFIGURATION FOR 
EXPERIMENTS 

Before proceeding to time varying currents, we note that 
the observations of Sec. II aid in designing other static con- 
figurations to produce a region of B=0, but A#0. 

Any localized current distribution Jj that produces a field 
Bj in a region where J^O can be transformed into a con- 
figuration that produces an A^O where B=0 by inventing a 
different current distribution that will produce a field equal to 
Jj. According to the hierarchy of Eq. (2.11), the needed cur- 
rent distribution is proportional to V X Jj . Since practical cur- 
rent sources are localized and produce a field outside them- 
selves, this shows that, by constructing the new current 
distribution a VxJ 1; there are many ways to produce a field- 
free region of space with A^O. 

The torus itself provides an example. The required current 
distribution is the curl of the torus current of Fig. 1, and is an 
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Fig. 3. Toms with counter-rotating surface current layers. The layer flowing 
in the -<p direction is on the outer surface of the tube; the layer flowing in 
the +<p direction is on the inner surface of the tube. 


azimuthally flowing double layer on the torus surface. The 
inner surface layer, say, flows in the direction <p , the outer 
surface layer in the direction -tp. B is in a thin sheet con- 
fined between the layers, and is in the toroidal direction, 
being a vector parallel to the current layer K of Fig. 1. As 
shown in Fig. 3 the vector potential circulates inside the 
torus, like the B of Fig. 1. This would provide a volume of 
slowly varying A. Unfortunately, this volume is physically 
inaccessible. 

However, if this torus is now cut through its limb at one 
place and then straightened out, we have a long narrow cyl- 
inder, of length 2 vb and radius a, with current flowing up its 
length on the outside cylinder surface, and back down its 
length on the inside. The magnetic field is circumferential 
about the axis, confined between the two current sheets. The 
vector potential is the same as the usual magnetic field of a 
solenoid. 

This configuration is itself topologically equivalent to the 
torus. Instead of changing the current direction and cutting 
the torus, we can stretch the torus of Fig. 1 in the z direction, 
its cross sectional circle of radius a being elongated into an 
ellipse, as in Fig. 4. After stretching to a length /, and taking 
a— >0, the resulting geometry is a “solenoid” of radius b and 
length /, with azimuthal magnetic field confined between 
the two axial surface current sheets. B=0 inside and outside 
this “solenoid.” Inside, at cylindrical radii p<b, A^O; out- 
side ( p>b ), A is only the “fringing field,” significantly dif- 
ferent from zero only near the ends. Inside the solenoid is a 
readily accessible volume for experiments. 

A inside is proportional to the flux of B, and so can be 
increased by thickening the walls, as in the cylindrical torus 
of Fig. 5. A inside, on or off axis, is axial and is, for /S2p, , 


A = 




(4.1) 



Fig. 4. Showing topological equivalence of torus and cylindrical “solenoid" 
by stretching the torus vertically. 



Fig. 5. Cylindrical torus with accessible region of large, constant A in the 
hole (at p<p/). 


A is maximized by increasing the ratio of outer to inner 
radius. For practical construction the logarithm is likely to be 
~1, providing a useable volume with radius p, of nearly 
constant A with magnitude as large as the largest occurring 
anywhere for any torus carrying the same current. If I varies 
in time, B always stays zero on axis, but E z =-SA z /dt is 
nonzero, providing a (smaller) region near the axis of van- 
ishing B, but nonzero A(/) and E(f). 

The value for A in (4.1) is, of course, gauge dependent, 
here being in the gauge V-A=0. As written this constant z 
component of A is V(zA), and so could be transformed 
away. But Eq. (3.8) requires that some A would then appear 
outside this cylindrical torus. 


V. QUASISTATIC FIELDS 


We return now to time dependent fields of a torus with 
circular cross section. 

If J and / vary in time, an observer in the near zone will 
see the previously static A vary as I(t)/r 3 and give rise to an 
electric field. If, further, VJ remains zero, the scalar poten- 
tial <f> vanishes, and the field is E=-<?A Idt. The quasistatic 
electric field pattern is the same as the static A, and, for 
harmonic variation with frequency w, is proportional to 
col/r 3 . A displacement current e^SE/dt— w 1 2 I/r 3 appears 
through, say, the area enclosed in the circular path drawn 
above the torus in Fig. 6. Then Ampere’s law 


1 8E 

VXB=p, 0 J+ ^2 


(5.1) 


requires there be a quasistatic azimuthal magnetic field to 
balance this displacement current. Integrating Eq. (5.1) over 
the area of radius p, the required magnetic field is given by 

1 d f 2 tt d f , 

2irpB=-^ — J dS-E--p-—J dp pE z . (5.2) 

Since E z ~ wl/p 3 for large p, this quasistatic B varies as 
co 2 I/r 2 for large r. 

This B must also be Vx A. But the just determined quasi- 
static B is not the curl of the quasistatic A, which vanishes. 
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VI. RADIATED FIELDS AND POTENTIAL 



Fig. 6. A quasistatic azimuthal magnetic field is required on the indicated 
loop due to the displacement current through it. 


The field discussion thus far therefore cannot be complete, 
and we need to look more closely at time-varying fields and 
potentials. 

In the Lorentz gauge, Maxwell’s equations reduce to the 
wave equation for the vector potential, whose Cartesian com- 
ponents are given exactly by the full retarded solution 


A(r,r) = 


M o 
477 



J(r',n 

Ir-r'l ’ 


where 


(5.3) 


|r-r' 

t' = t 

c 


(5.4) 


is retarded time. When A is expanded in powers of b/r, and 
each coefficient further developed in the low frequency ex- 
pansion in powers of kb=(ai/c)b, the lowest-order surviv- 
ing terms are given in Appendix A, Eq. (A17). The quasis- 
tatic vector potential is 

Aq S = e ~ V* f As tat ( r) , (5.5) 

where A^ is the static vector potential of Eq. (3.5). 

The quasistatic fields (lowest order in frequency) are 

fj, 0 VI ike 

Eqs=-; — j- (2 cos Or + sin 06), 

y 4 77 477 r 5 


Mo VI k 2 


(5.6) 


also proportional to torus volume. The quasistatic B does 
balance the displacement current as required by Ampere’s 
law, <?E QS /<?r = c"VxB 0S . As seen in the Appendix, E QS 
arises from the time derivative of the quasistatic potential 
Aq S , but B qs arises as the curl of the “inductive” A which is 
—ikr times Aq S . However, these quasistatic fields do not 
separately satisfy Faraday’s law, OB QS /Ot # - V X E qs . 
Rather, VxE QS =0, and <?B QS /dt is balanced by the curl of 
the part of E arising from the inductive A, which is kr times 
smaller than E QS . 

These fields, valid for kr< 1, are ordered according to 


VI k 2 

b qs p kr 


£qs<^£qs 
c c 


(5.7) 


For most charge-current distributions p, J, the quasistatic 
fields are simply the static fields with p(f), and J(t) in place 
of their static counterparts. This is decidedly not the case for 
a torus. 


Before computing the radiation fields from a torus we 
make some general observations. 

A. General comments on radiation 


In any gauge, the scalar potential is not needed for the 
radiated fields, since B=VxA, and E rad =-cnxB rad , where 
n is the outgoing unit vector. In the Lorentz gauge, A is 
given by (5.3). The radiated vector potential is the part that 
falls as 1/r, obtained by expanding l/|r — r'| and keeping only 
the 1/r term, 


A rad 


_M o_ 
477r 



J(r',/'). 


(6.1) 


If the time argument of J were t, the integral of J over all 
space at one instant of time would vanish when V-J=0. 
Retarded time means that wavelets emanating from different 
parts of the source do not quite cancel, allowing radiation. 
Retardation, of course, is the physical reason any divergence- 
less time-dependent current distribution of finite extent radi- 
ates. 

B rad is the 1/r part of VxA rad , or 
Brad^VxJrfV J(r\f'). (6.2) 


This does not vanish even for a torus. The integral depends 
on r only through t'. 


S3 


Vxj(r ',t') = Vt' x — = 

' dt' c |r-r 


1 r-r' OJ 

7T X ST. (6.3) 


where we have used the gradient of (5.4). Therefore, 


Brad 


_Mo_ 

477cr 




d\' 


r-r 

I r-r' 


03{r',t') 

dt' 


(6.4) 


The first factor in this integrand is (r-r')/|r- r'| 
-n -O(r'lr), so that 


Brad - ' 


Mo 

477cr 


nx 




(6.5) 


Now noting dt'/Ot = 1, the derivative may be pulled out, 


B, 


'rad“ 


Mo 

47rcr 


n x 4; d 3 r'5(r',t') = 

dt 


1 

— n x 
c 


^A rad 
dt ' 


(6.6) 


These equations also show E rad =-<?A rad A?r as it should. 

In any gauge, if A rad is nonzero, B rad and E rad are also 
nonzero. Unlike the static A, the radiated vector potential 
cannot, of course, be separated from its fields. That is, there 
is no such thing as a “radiated curl-free vector potential” 
occasionally referred to. It is easily shown that V-A rad =0, so 
that A,. ad is a transverse vector. 

The discussion so far is completely general, applying to 
any current distribution. 

As a rule, if a current increases in time from zero to a 
nonzero steady value, radiation is produced that leaves be- 
hind the static field of the nonzero current. When I was 
turned on for the torus, radiated A, E, and B propagated out, 
leaving behind a nonzero static A, and a static B which hap- 
pens to have value 0. The static A outside a torus is only 
accidentally curl-free because of high-symmetry geometry. 
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R 



divergence of the energy flux S=ExH (Poynting vector) 
enters the energy conservation law that follows from Max- 
well’s equations; VxS is arbitrary. S itself is arbitrary to the 
extent that the curl of any vector may be added to it. Yet we 
are quite accustomed to thinking of S as physically real. An 
arbitrary additive curl or gradient need not prevent a vector 
from being considered real. Konopinski 16 discusses how in 
classical electromagnetism A may be considered the poten- 
tial momentum per unit charge of a particle in external fields, 
just as <p is the potential energy per unit charge. 


Fig. 7. Resistive wire loop through a torus. 


As shown in Sec. VII it is more fundamentally to be consid- 
ered a transverse vector, with the defining property of being 
divergence-free everywhere. 


B. Physical argument for radiation 

It is instructive to understand physically why a system that 
encloses magnetic flux in an area that can be “looped” must 
radiate. Consider a resistive test wire that encircles the torus 
limb in a closed loop as in Fig. 7. Let R be the full resistance 
of the wire. When the torus current I varies, so does the 
magnetic flux <I> in the torus interior, and an emf 

W= <b d 1E=- — (6.7) 

J wire v* 

is produced, according to Faraday’s law. A current &IR flows 
in the test wire. The electrons in the wire know to move 
because the induced electric field drives them. 

All electromagnetic disturbances start where Vxj¥=0, in 
this case on the torus windings. The only way E can get to 
the wire is by propagating from the torus to the wire. This 
propagation, carried to larger distances, is radiation. 

If the torus and test wire are of very great radius L , then a 
propagation time of order Lie is necessary before E gets to 
the wire. But since Eq. (6.7) holds instantaneously, the radi- 
ated magnetic flux through the big wire loop is equal and 
opposite to the quasistatic flux inside the torus limb until the 
radiation front gets to the wire. It is only after the radiated 
fields pass that <J> in Eq. (6.7) can be taken to be the usual 
interior quasistatic flux. Only 1/r fields can account for the 
necessary emf. Radiation from a torus has been discussed by 
Baum. 14 

It is interesting to note that, from Eq. (6.7), and from Eq. 
(3.8), which holds for time varying conditions as well, the 
wire current is 


1 

R ~R 



1 

R 



( 6 . 8 ) 


so the total charge moved through the wire is 
f &{t) If 

(3(0= J dt= - ~ fdl-A. (6.9) 

That is, displaced charge is as much a direct measure of 
fdl-A as current is of <j>dl-E. From this point of view A 
could be considered just as real a physical variable as E. 
Feynman 15 stresses how A may be considered physically real 
in spite of the arbitrariness in its divergence. 

Parallel to the fact that only V xA enters Maxwell’s equa- 
tions and V A is arbitrary, it is worth noting that only the 


C. Explicit radiated fields 

To compute the radiated fields, we directly evaluate Eq. 
(6.1) for a harmonic toroidal current. For frequency w, J is 

/ . , 

J(r,0)= - — 8(s-a)e- ,a> ‘ a, (6.10) 

lirp 

where a is the unit vector in the direction of increasing a 
(Fig. 1). In Eq. (6.1) the s integral in d 3 r' = p' d<p' ds sda 
is trivial, leaving 

Un Id f 2 ™ f 2lr t 

A rad (r',t)= — - — dip' \ da a e~“°‘ . (6.11) 

47 t Inr Jo Jo 

In t ’ , |i — r' | is expanded 

r 1 . 

t'-t b-n-r '+0{r' 2 /cr) (6.12) 

c c 

so the integral in Eq. (6.11) becomes 

e~ io>T j dip' j da ae~' k r ', (6.13) 

where r=t — r/c is the observer’s retarded time, and k=(a>/ 
c)n. We have 

a= — (x cos ip' + y sin y/)sin a + z cos a (6.14) 

and 

k-r ' = k[(b + a cos a)cos(<p- <p')sin 0+a sin a cos ff], 

(6.15) 

where 9,ip are those of the observer. 

Equations (6.14) and (6.15) are to be inserted in Eq. 
(6.13). Rather than grapple with the resultant integrals we 
proceed with a low-frequency approximation, kb<& 1. When 
the exponential in the integrand in Eq. (6.13) is expanded, 
exp(-ik-r')=l-ik-r'-(k-r') 2 /2±..., it is easy to show the 
first two terms do not contribute to Eq. (6.13), and the re- 
maining integral is trivial. One finds, 

r 2 tt 

dip' 

Jo 


j- 2n 
JO 


da ae lk r 


= TT 2 k 2 ab sin 9(x cos 6 cos ip+y cos 9 sin ip 


— z sin 6) = TT 2 k 2 ab sin 99. 


(6.16) 


Then, if we call the time-dependent current / AC , to lowest 
nonvanishing order in kb the radiated A of a torus is, except 
for phase, 


A ra d 


k 2 I AC V 

47 r 47rr 


sin 99, 


(6.17) 


a factor of order (/ AC //)(A:r) 2 times the static A equation 
(3.5) in the same (Lorentz) gauge. Fields are 
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Fig. 8. Exaggerated model of vertical torus currents. 


, = Mo 

'rad a 
47 T 


Ic 3 IacV 

477 r 


sin d!p. 


®rad 


Mo ck 3 I AC V 

477 477r 


sin 00= — cr XB 


rad • 


(6.18) 


Radiated field amplitude is proportional to torus volume. 
The field pattern is that of an electric dipole. However the 
fields are proportional to a? rather than o> 2 . Physically, the 
origin of the fields is as follows. The torus of Fig. 1 has a 
current / flowing in the +z direction at p=b+a, and I flow- 
ing in the — z direction at p=b — a, with connecting paths 
above and below. Exaggerating the radii, Fig. 8 sketches a 
ring of upward flowing current balanced by a downward 
flowing current on axis. If I varies in time, both the center 
and the outer ring currents will radiate like individual electric 
dipoles. The electric dipole moment of the entire configura- 
tion is zero, since the charge displaced by one current is 
taken by the other. Since the sources are of different dimen- 
sions, their radiated fields do not cancel. Rather, the net field 
is proportional to ka, the separation relative to a wavelength. 
This explains the extra power of w and of a in the radiated 
fields relative to an ordinary electric dipole. 

One is reminded that even two simple dipoles, separated 
by d as in Fig. 9, together having zero dipole moment, will 
radiate an electric dipole field at frequencies a >^cc/d. At 
lower frequencies the quadrupole field dominates, the elec- 
tric dipole part being smaller by a factor of order wd/c. 


D. Multipole coefficients 

It is curious that all electric and magnetic multipole mo- 
ments of a toroidal current distribution vanish, and yet qua- 
sistatic fields and radiation do not. 

The general theory of multipole radiation [e.g., Jackson; 17 
Chap. 16] relates radiated fields to sources. The relevant 
source parameters are the electric and magnetic multipole 
coefficients, different from multipole moments in that they 
account for retardation across the source. Multipole moments 
form a complete set for static charge-current distributions, 
multipole coefficients for the time varying case. 

All magnetic multipole coefficients vanish for the torus, 
but the electric coefficients are nonzero. The electric coeffi- 
cients a E (l,m) are given in Jackson, 17 Eq. (16.91), and con- 
tain a contribution from charge density p, current density J, 
and intrinsic magnetization M. The term in p corresponds to 



Fig. 9. Two nearby small dipoles. 


the usual electric multipole moment when kr< 1 [generalized 
to include the radial radiation function ji(kr)], and vanishes 
for a torus, as does the term in M. 

The term in J is not zero. The extra factor kr in this term 
accounts for the destructive interference from opposite sides 
of closed loop currents. The lowest coefficient is Z=l, m= 0, 
and corresponds precisely to the previous explanation in 
terms of equal up and down currents at different radii. It is 
readily evaluated for the case kb< 1, and results in a radiated 
field equal to our Eq. (6.18). 

The general expression for a E (l,m) shows that a system 
with p=J=0, but with nonvanishing time varying M, also 
produces electric dipole radiation. 

There is no general formalism expressing quasistatic fields 
in terms of source moments as there is for static and radia- 
tion fields. By explicit calculation one finds quasistatic fields 
of a torus are nonzero and given by Eqs. (5.6). 


VII. GAUGE TRANSFORMATIONS AND SOURCES 

Due to the arbitrariness arising from a gauge transforma- 
tion, only certain parts of A are in principle measurable. One 
needs to inquire how A is related to current sources, since 
classically the sole raison d’etre of A is to have its curl equal 
to B. The longitudinal part of A can be combined with the 
scalar potential <f>, while the transverse part retains its iden- 
tity and is measurable. 

The source-free Maxwell equations (Faraday’s law, and 
the absence of magnetic monopoles) serve to define the fields 
in terms of the potentials, and the two equations with sources 
(Ampere’s and Gauss’ laws) become dynamical equations 
determining the potentials. If we make the transformation 

Ai-»A 2 =A 1 + V^, 

4>i~* <t>2= 4>\~ x, (7.1) 

where x > s an arbitrary function, gauge invariance assures 
that x drops out of the dynamical equations altogether and 
has no physical meaning. 


A. Transverse and longitudinal fields 

The relation VxA=B and the arbitrariness in V-A sug- 
gest writing all fields in their transverse (solenoidal) and lon- 
gitudinal (irrotational) parts. 18 Thus, for A, 

A= Aj-+ A l , 

VxA r =VxA=B, VxA t =0, (7.2) 

V-At-= 0, V-A t = V-A= arbitrary. 

If two vectors are equal, their transverse and longitudinal 
parts are separately equal. Then, decomposing E, B, and J 
into their parts, Maxwell’s equations can be written 

<?B r 

B t = 0, VXE 7 +— = 0; 


V-E L = p/e 0 , 

1 dY^j 

VxBr_ ?^r =/i ° Jr ’ 

with charge conservation 


(7.3) 

1 dE L 

~7-~dT = ^ L ’ 


dp 

V.J i+ - = 0. 


(7.4) 
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B. Maxwell’s equations without gauge transformations 


Now introduce the usual potentials A and <f>: 
B j -- V xAj- \ 


Ey' — 


dAj 

~dT’ 


E, 


dA L 

~dt 


V 4>. 


(7.5) 


Only A l and <p are affected by a gauge transformation, with 
Ejr left unchanged. A L can always be written in terms of a 
scalar ;y 0 , A L =V\o, so that E ; can be expressed as 

E L =-V^, (7.6) 

where 


A l , and the field associated with it, is a gauge artifice that is 
eliminated in Eq. (7.8), and has no physical meaning. 

The static A of a torus, Eq. (3.5), is related to the source 
current. This A is a transverse vector for which V -A=0 ev- 
erywhere and which happens to have VxA=0 outside the 
torus (but V xA=£0 inside). It is locally curl-free due to high 
symmetry. 

These equations make it clear that it is only the transverse 
part of the current source that produces radiation. This same 
result is apparent from the usual wave equation for B 

DB= pqV /r 0 V X J r (7.9) 

that follows from Maxwell’s equations. 


<A=<£+Xo (7.7) 

is a gauge invariant scalar potential with only an additive 
constant arbitrary. By itself A L is physically irrelevant, being 
a matter of gauge choice and having no relation to sources. It 
can be fully subsumed in the definition of the scalar potential 
in a gauge invariant way, and only <// and A r survive. In 
terms of them Maxwell’s equations are 

B r = V xA r , By = 0; 

ft A t 

E t = Ey=— V tjr, 



U\A r —-p () 3 T , — V — , (7.8) 

lit 6 q 

where n = V 2 -(l /c 2 )d 2 /dt 2 . 

Since A L has been combined with (p , the last equation for 
(// is redundant and may be dropped. It is equivalent to the 
remaining Poisson equation and the continuity equation 
(7.4). The Poisson equation for »/>, and wave equation for A r , 
are the same as the usual equations for (p and A in the Cou- 
lomb gauge, with the significant difference that now there 
has been no specification of gauge. Every quantity in Eq. 
(7.8) is gauge independent. 

With Maxwell’s equations in this form there is no longer 
any arbitrariness in the potentials, and the questions of gauge 
transformations or gauge invariance do not arise. Further, 
longitudinal sources, fields, and potentials completely de- 
couple from transverse ones. 

Maxwell’s equations are usually not formulated this way. 
Most theories insist on, and build in, relativistic invariance, 
and separately inquire as to gauge invariance. In contrast, the 
above equations build in gauge invariance (indeed, gauge 
irrelevance), but are not manifestly Lorentz invariant. The 
very division into transverse and longitudinal parts is not 
Lorentz invariant. However the formulation can be conve- 
nient for analysis of a given experiment in one reference 
frame. 

J £ gives rise only to if/. } T gives rise only to A T . Closed 
current loops, for example, which are commonly well ap- 
proximated as divergenceless, especially at low frequencies, 
are transverse currents, and produce only A r and transverse 
fields. In this same approximation, the complete fields of a 
torus with time-dependent current are purely transverse. 

The discussion following Eq. (6.6), that A^ cannot be 
separated from its fields, left open the possibility that the 
longitudinal part of A might be separable from its fields. But 


Vin. CONDUCTORS AS SHIELDS 

We inquire as to the effectiveness of conducting enclo- 
sures in shielding the vector potential. 

A. Boundary conditions on A 

Boundary conditions on A are conventionally derived 19 by 
applying the usual Gaussian pill box or Stokes loop over an 
interface between two media. But they can be obtained more 
directly from those on B by using the symmetry observed in 
Sec. II. 

Since the normal component B ± of B is continuous at an 
interface, and since this B is the A of a different problem, the 
normal component A ± of A is continuous also. Actually, for 
the general case, the discontinuity in A ± is formally gauge- 
dependent, since V - A need not vanish at the surface. But so 
long as the gauge is such that the volume integral of V-A 
across the surface remains zero, then A x is continuous. 

Similarly, that the tangential component B| is discontinu- 
ous by the normal integral of J means the tangential compo- 
nent A|| is continuous unless there is a delta function sheet of 
B on the surface. Excluding this unphysical condition, we 
have that all three components of A are continuous at an 
interface between two media with different e, p, and cr, even 
across a dielectric-conductor interface on which there may 
be a single-layer current sheet. 

The normal derivative of Aj| is not continuous if there is a 
skin current, and is determined by the usual boundary con- 
ditions on B|| . The discontinuity of E L at a surface charge 
density shows up in the potentials as a discontinuity in the 
normal component of Vcp. 

B. Conductor moving through field-free region 

If a conductor is moved through a region of space where 
E=B=0, but A^O, there are no physical effects, for there 
are no forces to produce any. Since A is continuous across 
the conductor surface, A effectively penetrates freely through 
the conductor. 

To an observer at rest on the conductor, A is changing in 
time at a rate dA]dt=\-VA, where v is the conductor veloc- 
ity. By a Lorentz transformation, the observer sees a scalar 
potential <p=—\-A, just such as to keep E=-<?A/<?r-V</> 
equal to zero. 

C. Shielding by a perfectly conducting enclosure 

Consider a localized current distribution with its attendant 
vector potential with B=0. We enclose it in a conducting 
shield. Two “experiments” are considered: 
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Shield 



Fig. 10. Torus shielded by a larger torus. 


(1) The current and vector potential are pre-established and 
the shield box is erected around the source in the space 
where A=£0. 

(2) The shield box is constructed while the source is off, 
then the source is turned on. 

We will choose a torus for the source. 

1. Pre-established dc field 

With a steady current flowing, A is nonzero outside. As 
the conductors are assembled, they move in the space 
through A. A passes freely through the metal pieces; there is 
no interaction between A and the conductors. We end up 
with the torus with its unperturbed static A, enclosed by a 
conducting box, with A^O inside and outside the box. A has 
no way of knowing that the shield box was built. 

2. dc field turned on after shield is built 

In this case we erect the shield box around the dead torus, 
then turn the current on. The shield interrupts A, E, and B 
that are propagating out. So that concepts are not clouded by 
noncompatible geometries, choose for the box a concentric 
torus. 

E, B, and A remain zero inside the shield metal. There- 
fore, on the inner surface of the shield, E, A, and B ± also 
remain zero. Surface currents which terminate By are gener- 
ated on the inner shield surface. The space between the driv- 
ing torus and the shield torus is filled with B and A. 

If / is suddenly turned on and then held constant, the final 
static configuration is as shown in Fig. 10. A skin current 
flows on the inner surface of the shield, in the opposite di- 
rection of the driving current. An azimuthal B field persists 
within the shield box, opposite the main B in the driving 
torus, so that the total flux through the large torus arm stays 
zero. The static A is in the toroidal direction between the two 
tori, but is 0 at the shield. A, E, and B all vanish outside the 
shield. The box shields A as well as the fields. 

There is thus a significant difference according as the 
shielding enclosure is constructed after or before the driving 
current is activated. In the former case, the static A is non- 
zero outside the shield, but in the latter case it stays zero. 

If the enclosing box is a finite conductor, currents and 
fields diffuse into the metal and penetrate the shield of thick- 
ness d after a time of order 

t d =/x 0 crd 2 (8.1) 



Fig. 11. Hierarchy of current distributions based on torus. 


(^*=180 /jus for 2 mm thick Al). After this time, A will ap- 
pear outside the “shield box,” and the exterior dc vector 
potential will have been established. 


IX. VARIATIONS ON A THEME 

The observations in Sec. II lead to a physically meaningful 
hierarchy of static current configurations. 


A. Hierarchy of J, B, A 

Consider the current distribution of the torus, Fig. 1, re- 
drawn in Fig. 11(2). Call its current J 2 , field B 2 , and vector 
potential A 2 . 

Construct a new current distribution J 3 =B 2 having only an 
azimuthal (<p) component. This is the current of an ordinary 
current loop (except that J 3 drops off as 1/p within the loop). 
Its field B 3 will be B 3 =A 2 . Its vector potential must be com- 
puted anew. A 3 is also azimuthal, and lines of A 3 form circles 
about the z axis in space, as sketched in Fig. 11(3). 

Continuing, construct J 4 =B 3 . This current distribution 
would be that of the discharge current of a battery immersed 
in a partially conducting fluid. The associated magnetic field 
is B 4 =A 3 , being circles about the axis, encircling lines of J 4 . 
The vector potential A 4 would have a field line pattern close 
to that of J 4 . 

Proceeding in the opposite direction, from configuration 
(2) invent a current distribution Ji<*VxJ 2 whose field 
Bj=J 2 , and vector potential A 1 =B 2 . This is the azimuthal 
double-layer current sheet discussed in Sec. IV, whose field 
and vector potential vanish outside the torus. It is sketched in 
Fig. 11(1). 

One could go one step further and construct four alternat- 
ing layers of toroidal currents, J 0 , confining B 0 to be the 
same double layer as Jj . Aq would be like Bj , confined to a 
thin toroidal sheet, the same as J 2 . A, B, and J vanish inside 
and outside the torus. 

Developments similar to this can be based on a common 
cylindrical solenoid or any other current distribution. The 
hierarchy is summarized by 
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( 10 . 1 ) 


# 



Fig. 12. Toroidal coax line. 


J m •••= Vxj t 
B m . . . = J ! 

A m ---~ ®1 

A, 


V xj 2 

J 2 = V x J 3 

B 2 = J3 = 

A 2 = B 3 — 

A 3 


which is written more succinctly in Eq. (2.11). 


...Vxj„ 

-J. 

•••A„ 


B. Toroidal coax 

In the configuration of Fig. 11(1), one can separate the 
inner current sheet from the outer one, collapsing it to an 
inner wire ring concentric with the torus limb, as in Fig. 12. 
One has a toroidal coax line. J and A are azimuthal, A being 
nonzero both inside the inner conductor, and between the 
inner and outer one. B is in the “toroidal” direction, confined 
to between the inner conductor and outer “shield,” being the 
magnetic field of an ordinary coax line. A vanishes outside 
the entire configuration and inside at the outer conductor. 

As the inner wire is an ordinary current loop, this configu- 
ration is a “shielded current loop.” The shield is an ordinary 
(perfect) conductor and shields both A and B. The same is 
true for ordinary straight coax cables. 


X. A SYMMETRY OF MAXWELL’S EQUATIONS 

It is interesting to inquire whether the magnetostatic ob- 
servations of Sec. II on the J, B, A hierarchy can be extended 
to time varying currents and fields. 

Due to the general relations B=VxA, V A=0, and the 
static relation J=VxB, as well as V-B=0, we have for the 
arbitrary time dependent case. 

The vector potential A 2 (r,f) of any current distribution 
J 2 (r,f) with field B 2 (r,f) is the same as the instantaneous 
static magnetic field B,(r,t) of a current J x (r,t) instanta- 
neously equal to B 2 (r,t). 

Here A 2 is in the Coulomb gauge. It appears this rule cannot 
be extended to develop a hierarchy of J, B, A configurations 
in parallel with the static case since the displacement current 
destroys the analogy. But a closely related hierarchy still 
holds. 

In the Lorentz gauge Maxwell’s equations are 


dA 

B= V xA, E= - — — V<£, 

at 


□A=— fi 0 J, □</>=-—. 

e o 

Let these hold for a current J, with fields Ej , , and poten- 

tials Aj , <fi] . 

Form the new current 

J 2 = VxJi, V-J 2 = 0. (10.2) 

Clearly, its potentials are 

A 2 = V x Aj = Bi , 

<f> 2 = 0. (10.3) 

The associated fields are 



~~dt 


= VxE 1; 


/ dEj) 

B 2 — V xA 2 — V xBj — /z 0 ( Ji + e 0 ~^~J ■ 

In summary, 


/ J \ 

/ J \ 

B 1 

B 1 


1 E 

\ A / 

\ A / 


(10.4) 


(10.5) 


a complete parallel to the magnetostatic case (2.11). Even for 
full time dependence the curl of one set of currents, fields, 
and potentials is another set, so that a solution for one J 
provides the solution for other problems. That this should be 
true for Maxwell’s equations is less obvious than for magne- 
tostatics, although it follows as well simply by taking the 
curl of Maxwell’s field equations. It differs from magneto- 
statics in that E^0, and VxB is no longer J, so the hierarchy 
is not as tight. This is made more explicit by writing Eq. 
(10.5) as 


/ h) 


/ Vxj,\ 

1 ° \ 

B 2 


J , 1 

*Ei Idt 

Ej , 

— 

1 0 

+ -dBi/dt 1 

w 


\ Bj 1 

\ 0 / 


( 10 . 6 ) 


in which we have dropped factors of and in the second 
line for simplicity. A problem-specific numerical coefficient 
with dimension of length should also appear multiplying the 
right-hand sides of Eqs. (10.5) and (10.6). Based on any 
given J(r,f) one can develop a hierarchy of physical current 
and field distributions in parallel with the magnetostatics 
case. 


A. Fields of a rotating torus 

One instance in which this symmetry is useful is in com- 
puting the fields of a rotating torus. 

Let a dc current flow in the torus windings of Fig. 1. 
Suppose the current is delivered to the windings by slip rings 
so the torus is free to rotate about, say, the y axis, maintain- 
ing its dc current. (Or we could imagine a superconducting 
arrangement, with no need to drive the currents once estab- 
lished.) When stationary, there is a vector potential, but no 
fields, outside. When the torus is spinning, are there external 
fields? 
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Clearly there are, for to an external observer the static 
A also rotates, producing nonzero E=-<2A jdt, without a 
canceling — V0. Since GE/dt is also nonzero, 
VxB=(l/c 2 )o!E/<?f#0, and therefore B does not vanish out- 
side. An oscillating EMF is generated around a fixed loop in, 
say, the jc=0 plane that encircles the torus limb in its original 
position. A rotating torus carrying a dc current possesses 
nonzero near zone fields. It also radiates. The question is 
how to calculate the fields. The symmetry discussed above is 
quite powerful in this regard, and makes the task easy. 

The surface current density on the torus is the curl of the 
current density of an ordinary current loop with the same 
dimensions. According to Eq. (10.5) then, the exact fields 
and potential of the rotating torus are precisely the curl of the 
corresponding fields and potential of a spinning current loop. 
These latter fields are not difficult to compute. Here we 
sketch only the radiation fields. 

More precisely, a torus wound with N turns carrying total 
current I=i JV, has a current density J proportional to the 
curl of the current density J' of a loop carrying current 
I' = Tra 2 J'\ 


J=kVxJ', 


a 2 I 
2b I 7 ' 


(10.7) 


The coefficient k should appear multiplying the right-hand 
side of (10.5) and (10.6). 

A loop of radius b, of thin wire with radius a<b, carrying 
current /' has a magnetic moment m = Trb 2 !' . When rotating 
with angular velocity <u in the positive sense about y its 
components are 

m z =m cos (at, m x — m sin (at. (10.8) 

Each of these is in turn the magnetic moment of a stationary 
loop with oscillating current. Each radiated electric field is of 
the standard form, and the total is the superposition of the 
two: 


Er s =KVxE‘°J 


/i 0 c Vk 3 I e i< ' kro>l> 

Att Att r 


[-sin 6 x 8 x +i sin 8 Z 8 Z ], 


(10.12) 

k has been replaced using Ktn = (ira 2 b 12)1 = (VIA tt )1 . The 
full radiation pattern can be mapped from Eq. (10.12). As for 
a loop, this is the field of two stationary perpendicular 
toruses carrying oscillating currents, each with fields of the 
form (6.18). 

For an observer on the rotation axis (+y), for example, 
sin 0 x =sin 8 Z = 1, 8 X = —x, 6 Z = -z, and 


ES“(x=0,y = r,z=0) = 


/m 0 c Vk 3 I e 
477 477 


i(kr — <ot) 


[jc — iz]. 
(10.13) 


Radiation in the +y direction is right-circularly polarized. 

One can similarly compute the quasistatic fields of a spin- 
ning torus from those of a spinning loop. 

The energy of the torus consists of its rotational kinetic 
energy plus the magnetic field energy inside; these supply the 
energy radiated. The kinetic energy is an artifact of the mass 
of material chosen for fabrication, and can in principle be 
made as small as desired. Therefore, the energy radiated 
comes from the enclosed magnetostatic field energy. Due to 
radiation the dc current of a spinning torus will decay. 

This example has been illustrative only. One can employ 
the symmetry noted in this section to problems less academic 
than a rotating torus. 
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ploop_p(jr) , p(z) 

^rad “ *\ad ^ ^rad > 

, , Un (a 2 m ... 

E$= -i ~ e i(kr-a>t) sin g * 

rad 477 cr x * x 


(10.9) 


fz 0 (a m 


E (z) = — 

trad 4 77 cr 


i(kr~(ot ) 


sin 0 z (p z , 


where the real part is understood. Here 8 Z and <p z are the 
usual polar and azimuthal angles in spherical coordinates 
defined on the z axis, and G x and (p x are those of a spherical 
system defined on the x axis. The magnetic field is 

B‘“P= ^ fxE^. (10.10) 

For these, as for any radiation fields, 

<?B rad u> 

V xE rad = — = i«B rad = i — rxE rad . (10.11) 

Then, using r X ip x = - 0 X , etc., the radiated field of a spin- 
ning torus comes out 


APPENDIX A. EXACT FIELDS OF A TORUS 


The complete vector potential and fields of a harmonically 
driven torus may be computed to arbitrary accuracy by ex- 
panding in small parameters. 

The full expression for A(r,t) in the Lorentz gauge is Eq. 
(5.3). Taking 


J(r' ,t’) = e- ,w '' J(r' ) = e -<^ e ‘*l r - r ' lj( r ' ) ( A l) 

and A=A a ,(r)e -l “', then 

r e **|r— r'l 

TTTT J ( r ')- ( A2 > 

} r— r 


A w (r) = 


Mo 

477 


V-A=0, and since V-J=0, the scalar potential vanishes. We 
consider only r>b and kb< 1, but expressions will be valid 
in the near zone &r<l, or far zone kr> 1. When kb^l, 
azimuthal asymmetries arise due to propagation time delays 
around the torus when driven at one point. We do not take 
into account these asymmetries. 

The dependence on r and r' is separated using the stan- 
dard expansion 


e (*|r-r'| “ 

1 77 = ^2 (2/+1) jAkr')h { ;\kr)PAfi), r>r', 

|r-r | /=o 

(A3) 


727 


Am. J. Phys., Vol. 63, No. 8, August 1995 


N. J. Carron 


727 



where is the spherical Bessel function, P/ is the Legendre 
polynomial, and /i=cos(r,r'). Powers of kr are explicit be- 
cause the Hankel function 


h ( }\kr) = 


1 e 


ikr S 


(S+q)! 


Y ikr q =o \2kr 


Then, 


A --E e ~ 2 0 SAt) £MY • (A5) 


ikr 


1 V 


where 


S/k)=f d 3 r'j 


jAkr')PAv)J(r') 

characterizes the source and observer angle, and 
(2/+1) (S+q)l 


f/q — 


Y~< 2 *q\(/-q)\ 


(A4) 


(A6) 


(A7) 


Since for small argument j Akr')={kr'Y /{ 2/+1)!!, 
Equation (A6) shows 


SAk) = (kbfG / , />0 


(A8) 


with G/ = G/(k) independent of k to lowest order. The only 
exception is S ^{kby. Using this in (A5) and interchanging 
summations gives 


Mo e ikr " 

A w =t-— /ooSo+X 

47T r I q = o 


, ( kbY 

G/// > oup 


(A9) 


where the prime means the /=0 term is deleted from the 
sum. Now setting A—q+j, 
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This is an explicit series in powers of b/r with coefficients 
that are rapidly converging series in kb. 

It is not difficult to show that the combination of azi- 
muthal symmetry and reflection symmetry in the z=0 plane 
implies 


S/{k) = GAk) = 0, / odd 
so that only even indices survive in (A10): 
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A w contains all powers of 1/r, but as >0 the static vector 
potential 
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contains only odd inverse powers of r. 
The first few terms of (A12) are 
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In evaluating the higher powers, the k dependence of G/ 
must be included. G 2 has already been evaluated, since ac- 
cording to (A13) it is determined by the static 1/r 3 potential, 
which was obtained in Sec. III. Comparing the first term of 
(A13) with (3.5) gives 
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and direct evaluation shows 
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Then, taking the /’ s from (A7), one gets 
Mo VI 
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plus smaller terms of higher order in b/r and kb. Here 
V=2n 2 a 2 b= toms volume. 

The first term (1) in square brackets is the lowest order 
static vector potential, separately obtained in Eq. (3.5). The 
ikr term is the “inductive” potential. The term in k 2 r 2 is the 
radiated potential, previously obtained in Eq. (6.17). 

The fields to the same order are 
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E— — — = i o>A, 
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(Ai8) 


having terms behaving as a )I/r 3 , w 2 I/r 2 , and w 3 I/r; and 
Mo VI . , , , k 2 

B= V XA= P 2 — e- ,a{, - T/c) — [1- ikr] sin 0q>, 

4 7T 4 7T r 

(A19) 

behaving as uTl/r 2 and a?I/r. 

In the near zone, the quasistatic magnetic and electric 
fields, given in (5.6), are related by Eq. (5.7). 
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Teaching physics with 670 nm diode lasers — construction of stabilized 
lasers and lithium cells 

K. G. Libbrecht, R. A. Boyd, P. A. Willems, T. L. Gustavson, a) and D. K. Kim b) 

Norman Bridge Laboratory of Physics, California Institute of Technology 12-33, Pasadena, 

California 91125 

(Received 23 September 1994; accepted 22 February 1995) 

We describe the construction and operation of stabilized 670 nm diode lasers for use in 
undergraduate teaching labs. Because they emit low-power visible radiation, 670 nm lasers are safe 
and aesthetically pleasing, and thus are an attractive alternative to near-infrared diode lasers in the 
undergraduate laboratory. We also describe the fabrication of a robust and reliable lithium atomic 
vapor cell, which can be used with the 670 nm diode lasers to perform a variety of atomic physics 
experiments. © 1995 American Association of Physics Teachers. 


I. INTRODUCTION 

As inexpensive sources of narrowband tunable coherent 
light, semiconductor diode lasers are becoming important re- 
search tools with widespread applications in modem aca- 
demic and industrial laboratories. Because of their relatively 
low cost (compared to other tunable laser sources) they are 
also well-suited for incorporation into undergraduate teach- 
ing labs. By scanning the laser emission over atomic or mo- 
lecular resonance lines, a variety of interesting and funda- 
mental physics experiments can be performed. 


A recent paper by MacAdam et al . 1 describes several tech- 
niques for constructing and operating stabilized diode lasers 
for the undergraduate laboratory. These lasers operate in the 
near infrared, at 780 or 852 nm, and can be tuned to scan 
over various rubidium or cesium resonance lines, respec- 
tively. Unfortunately, radiation at these wavelengths is al- 
most completely invisible to the human eye, which creates 
several problems when incorporating near-IR lasers into un- 
dergraduate labs. A primary concern is safety, since even a 
milliwatt of laser radiation can cause irreparable eye damage. 


729 


Am. J. Phys. 63 (8), August 1995 


© 1995 American Association of Physics Teachers 


729 



